Astronomy & Astrophysics manuscript no. Khelifi 


February 2, 2(?08 


(DOI: will be inserted by hand later) 





Discovery of extended VHE gamma-ray emission from the 
asymmetric pulsar wind nebula in MSH 15-52 with H.E.S.S. 

F. Aharonian 1 , A.G. Akhperjanian 2 , K.-M. Aye 3 , A.R. Bazer-Bachi 4 , M. Beilicke 5 , W. Benbow 1 , 

D. Berge 1 , P. Berghaus 6 *, K. Bernlohr 1,7 , C. Boisson 8 , O. Bolz 1 , I. Braun 1 , F. Breitling 7 , A.M. Brown 3 , 
iy) ■ J. Bussons Gordo 9 , P.M. Chadwick 3 , L.-M. Chounet 10 , R. Cornils 5 , L. Costamante 1,20 , B. Degrange 10 , 

A. Djannati-Atai 6 , L.O'C. Drury 11 , G. Dubus 10 , D. Emmanoulopoulos 12 , P. Espigat 6 , F. Feinstein 9 , 
^ ! P. Fleury 10 , G. Fontaine 10 , Y. Fuchs 13 , S. Funk 1 , Y.A. Gallant 9 , B. Giebels 10 , S. Gillessen 1 , 

J.F. Glicenstein 14 , P. Goret 14 , C. Hadjichristidis 3 , M. Hauser 12 , G. Heinzelmann 5 , G. Henri 13 , 
Q_i. G. Hermann 1 , J. A. Hinton 1 , W. Hofmann 1 , M. Holleran 15 , D. Horns 1 , O.C. de Jager 15 , B. Khelifi 1 , 
Nu. Komin 7 , A. Konopelko 1,7 , I.J. Latham 3 , R. Le Gallou 3 , A. Lemiere 6 , M. Lemoine-Goumard 10 , 
N. Leroy 10 , T. Lohse 7 , O. Martineau-Huynh 16 , A. Marcowith 4 , C. Masterson 1 ' 20 , T.J.L. McComb 3 , 
M. de Naurois 16 , S.J. Nolan 3 , A. Noutsos 3 , K.J. Orford 3 , J.L. Osborne 3 , M. Ouchrif 16 ' 20 , M. Panter 1 , 

G. Pelletier 13 , S. Pita 6 , G. Piihlhofer 1 ' 12 , M. Punch 6 , B.C. Raubenheimer 15 , M. Raue 5 , J. Raux 16 , 
S.M. Rayner 3 , I. Redondo 10 ' 20 **, A. Reimer 17 , O. Reimer 17 , J. Ripken 5 , L. Rob 18 , L. Rolland 16 , 

G. Rowell 1 , V. Sahakian 2 , L. Sauge 13 , S. Schlenker 7 , R. Schlickeiser 17 , C. Schuster 17 , U. Schwanke 7 , 
M. Siewert 17 , H. Sol 8 , R. Steenkamp 19 , C. Stegmann 7 , J.-P. Tavernet 16 , R. Terrier 6 , C.G. Theoret 6 , 
M. Tluczykont 10 ' 20 , G. Vasileiadis 9 , C. Venter 15 , P. Vincent 16 , H.J. Volk 1 , and S.J. Wagner 12 



in 



2 
3 



o 



CZ5 



O 

in . 

Max-Planck-Institut fiir Kernphysik, Heidelberg, Germany 
Yerevan Physics Institute, Armenia 
University of Durham, Department of Physics, U.K. 
C^h' 4 Centre d'Etude Spatiale des Rayonnements, CNRS/UPS, Toulouse, France 

Universitat Hamburg, Institut fiir Experimentalphysik, Germany 
APC, Paris, France *** 

Institut fiir Physik, Humboldt-Universitat zu Berlin, Germany 
8 LUTH, UMR 8102 du CNRS, Observatoire de Paris, Section de Meudon, France 

Groupe d'Astroparticules de Montpellier, IN2P3/CNRS, Universite Montpellier II, France 
Laboratoire Leprince-Ringuet, IN2P3/CNRS, Ecole Polytechnique, Palaiseau, France 
, 11 Dublin Institute for Advanced Studies, Ireland 

5^ ' 12 Landessternwarte, Konigstuhl, Heidelberg, Germany 

13 Laboratoire dAstrophysique de Grenoble, INSU/CNRS, Universite Joseph Fourier, France 

14 DAPNIA/DSM/CEA, CE Saclay, Gif-sur-Yvette, France 

15 Unit for Space Physics, North- West University, Potchefstroom, South Africa 

16 Laboratoire de Physique Nucleaire et de Hautes Energies, IN2P3/CNRS, Universites Paris VI & VII, France 

17 Institut fiir Theoretische Physik, Lehrstuhl IV: Weltraum und Astrophysik, Ruhr-Universitat Bochum 
Germany 

Institute of Particle and Nuclear Physics, Charles University, Prague, Czech Republic 



9 
10 



18 



19 University of Namibia, Windhoek, Namibia 



European Associated Laboratory for Gamma- Ray Astronomy, jointly supported by CNRS and MPG 



Accepted by A&A 



Abstract. 

The Supernova Remnant MSH 15-5,2 has been observed in very high energy (VHE) 7-rays using the H.E.S.S. 
4-telescope array located in Namibia. A 7-ray signal is detected at the 25 sigma level during an exposure of 22.1 
hours live time. The image reveals an elliptically shaped emission region around the pulsar PSRB1509-58, with 
semi- major axis ~6' in the NW-SE direction and semi-minor axis ~2'. This morphology coincides with the diffuse 
pulsar wind nebula as observed at X-ray energies by ROSAT. The overall energy spectrum from 280 GeV up to 
40TeV can be fitted by a power law with photon index r=2.27±0.03 s tat±0.20 S yst- The detected emission can be 
plausibly explained by inverse Compton scattering of accelerated relativistic electrons with soft photons. 



Key words, gamma-rays: observations - ISM: supernova remnants - pulsar wind nebula - ISM:individual objects: 
PSRB1509-58 - ISM:individual objects:G 320.4-1.2 



1. Introduction 

The Supernova Remnant (SNR) MSH 15-52 (G 320.4- 
1.2), first observed as an extended non-thermal radio 
source by Mills et al. l(llT61 ). is a complex object with 
an unusual morphology. Radio observations by Caswell et 
al. I|1981[) reveal a roughly circular SNR ~30' in diam- 
eter with a bright feature in the NW rim and a fainter 
one in the SE. The 10' diameter NW source (G 320.4-1.0) 
coincides with the H a nebula RCW89. 

Einstein X-ray observations of MSH 15-52 by Seward 
and Harnden (1982) led to the discovery of the 150- 
ms pulsar PSRB1509-58 located within the SNR shell, 
surrounded by a diffuse extended non-thermal compo- 
nent. The existence of a pulsar wind nebula (PWN) was 
later confirmed by ROSAT IjTrussoni et al. 19961 Brazier 
& Becker UnnZJ and ASCA (|Tamura fit al. 1996) as an 
extended emission region surrounding the pulsar with 
a power-law photon index of ~2.0. The PWN mor- 
phology is clearly visible in the ROSAT PSPC data 
(Trussoni et al. 1996) which show an elongated structure 
roughly centered on the pulsar with two arms extend- 
ing several arcminutes along the NW and SE directions. 
These features were more recently confirmed by detailed 
Chandra observations (Gaensler et al. [2002 ) . 

Gaensler et al. 1|1999|) . comparing X-ray and radio ob- 
servations, concluded that the radio SNR MSH 15-51? and 
the pulsar PSRB1509-58 are parts of a single system in- 
teracting via a pair of opposed collimated outflows (i.e. 
the PWN). The distance to the object and the pulsar spin- 
down age were estimated to 5.2±1.4kpc and 1700 years, 
respectively. 

The PWN hard X-ray spectrum was measured by 
BeppoSAX ( Min eo et al. 200T)l and is fitted by a power 
law with photon index of 2.08±0.01 in the 1-200 keV 
range. Interpreted as synchrotron emission from electrons 
within the PWN, these results suggest the presence of 
electrons accelerated up to energies of tens of TeV, pos- 
sibly leading to VHE 7-ray emission through an inverse 
Compton (IC) process as suggested by du Plessis et al. 
I|1995[) . Early observations by the CANGAROO experi- 
ment (jSako et al. 2000|) yielded a marginal 7-ray signal 
above 1.9 TeV corresponding to ~10% of the flux of the 
Crab Nebula. 

2. H.E.S.S. Observations and Results 

MSH 15-52 was observed using the atmospheric 
Cherenkov technique with the full 4-telescope H.E.S.S. 
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(High Energy Stereoscopic System) array (jHinton 2004|l . 
The observations were made between March and June 
2004 in wobble mode, for which the target position is off- 
set by 0.5° with respect to the centre of the field of view. 
Only data with good atmospheric conditions were selected 
yielding a total live time of 22.1 hours. The selected data 
were taken at a mean zenith angle of 37°. A standard 
analysis is used for the selection of 7-ray candidates, for 
the stereoscopic direction reconstruction of showers, and 
for spectral analysis iAharoni an et al. 2005J) . resulting in 
an average energy threshold of ^280 GeV at this zenith 
angle. For each sky bin (0.6'x0.6'), the background level 
is estimated from events falling in a ring centered on this 
test position with a mean radius of 1° and with an area 7 
times that of the cell. 

An excess with a significance of ~25 a (using the 
likelihood method of Li & Ma I1983[) is detected within 
the region of radius 0.14° centered at the pulsar posi- 
tion with the standard point source analysis. Figure ^ 
shows the excess map of 7-ray candidates with image 
sizes greater than 400 photoelectrons (p.e.). The latter 
cut is meant to reduce the angular resolution to <0.07° 
(70% containment radius) and raises the energy threshold 
to ~900 GeV. The map is convolved with a Gaussian of 
(7=0.04° in order to smooth out statistical fluctuations. 
An extended emission is seen along the jet axis of the pul- 
sar UTrussoni et al. 19~96)l in the NW-SE direction. There 
is no evidence for VHE 7-ray emission from the shell of 
SNR G 320.4-1.2, including the NW part (RCW89) from 
which a combination of thermal and non-thermal emis- 
sion was seen in X-rays by ASCA IjTamura et al. 1996) 
and Chandra (Gaensler et al. 2 002fl . 

The unsmoothed excess map (with image sizes above 
400 p.e.) is fitted by a two-dimensional Gaussian con- 
volved with the point spread function (PSF) in order to 
extract the intrinsic dimensions. The best fit position of 
the Gaussian centroid is (15hl4m7s±21s, -59°9'27"±11") 
which is displaced by more than 3 a from the pulsar po- 
sition (15hl3m55.6s, — 59°8'8.9") taking into account the 
20" systematic error on the pointing in each direction. The 
fitted direction of the major axis is 41°±13° with respect 
to the RA axis, and is compatible with the X-ray direction 
(60° ±5°) observed by Chandra. The intrinsic standard 
deviations along the major and minor axis are 6.4'±0.7' 
and 2.3'±0.5', respectively. These dimensions are compa- 
rable to t hose of the diffuse PW N as observed by ROSAT, 
~10'x6' l|Trussoni et al. 1996|) . Figure |21 shows the pro- 
jection along the major and minor axis of the count map 
relative to the pulsar position and illustrates the morpho- 
logical features described above. The transverse profile is 
approximately symmetric. The longitudinal profile, in con- 
trast, appears to extend more in the SE than in the NW 
direction relative to the pulsar position. Given the spatial 
coincidence and the similar morphology of the H.E.S.S. 
picture and the X-ray data, we identify this new VHE 
7-ray source with the pulsar wind nebula in MSH 15-52. 

The total 7-ray excess within a circle centered on the 
Gaussian centroid with a radius of 0.3°, which encom- 
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Fig. 1. Smoothed excess map from MSH 15-5,2 in arbitrary units (a.u.). The map is smoothed with a Gaussian of 
(7=0.04° and only events with image sizes above 400 p.e. are used in order to improve the H.E.S.S. angular resolution. 
The white contour lines denote the X-ray (0.6-2.1 keV) count rate measured by ROSAT (|Trussoni et al. 1 996). The 
black point and black star lie at the pulsar position and at the excess centroid, respectively. The right-bottom inset 
shows the simulated PSF smoothed identically. 

passes the entire VHE 7-ray emission, is 3481±129 events 
after standard cuts. Fig. shows the reconstructed spec- 
trum from these events. The data are consistent with a 
power law (dN/dE(xE~ r ) up to ~-40TeV with a photon 
index of r=2.27±0.03 sta t±0.20 syst and a differential flux 
at 1 TeV of (5.7±0.2 stat ±1.4 syst ) x 10~ 12 TeV^ 1 cm" 2 s" 1 . 
The x 2 /d.o.f. of the fit is 13.3/12. The corresponding in- 
tegral flux above 280 GeV represents ~15% of the Crab 
Nebula flux above the same threshold. The flux distribu- 
tion, as measured on a run by run basis, is compatible 
with steady emission (the x 2 /d.o.f. is 47.0/51). 



3. Discussion and Conclusions 

The H.E.S.S. observations of MSH 15-52 provide the first 
image of an extended PWN in the VHE range. The emis- 
sion region is clearly extended along the jet axis. If this 
VHE 7-ray emission arises from inverse Compton (IC) 
scattering of target photons by electrons, the latter should 
emit synchrotron radiation below 1 keV. Thus, the mor- 
phology similar to that seen in soft X-rays by ROSAT sug- 
gests a leptonic origin of the H.E.S.S. signal. Assuming 
the target photons are uniformly distributed, IC emis- 
sion directly reflects the distribution of high-energy elec- 
trons, unlike synchrotron emission which also reflects spa- 
tial variations of the magnetic field. 




Fig. 2. Projection of the unsmoothed excess map along 
the major (top) and minor (bottom) axes relative to 
the pulsar position. The dot-dashed lines are the best- 
fit Gaussians to these distributions, and the dashed lines 
show the distribution for a Monte-Carlo point source ex- 
cess at the pulsar position. 

The lack of VHE 7-ray emission in the NE-SW direc- 
tion shows that the elongated morphology of the PWN 
in X-rays is not due to a lower magnetic field strength, 
B, further away from the jet axis. However the reason for 
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Fig. 3. Reconstructed VHE 7-ray spectrum of the whole 
nebula. The lines are the best model fit of the one zone 
model described in section 01 

the apparent lack of high-energy electrons in the equato- 
rial direction is not clear. If the spin axis is aligned with 
the magnetic axis of the pulsar, as suggested by Brazier 
& Becker l)1997|) . Sulkanen & Lovelace (JT990) show that 
most of the energy, angular momentum and electric cur- 
rent are carried within the jet, which may enhance the 
accelerated electron density wihtin jets. The VHE 7-ray 
emission asymmetry along the jet axis relative to the pul- 
sar may be interpreted as Doppler boosting within the 
two opposite jets < |Pelling et al. 1987| ). The two jets may 
also be confined differently by the surrouding medium, as 
suggested by the appearance of RCW 89 to the NW. 

The total power, L, radiated by MSH 15-52 in the en- 
ergy band 0.3-40 TeV is ~1.0 x 10 35 ergs -1 (here we adopt 
a distance of 5 kpc) compared to ~8.0 x 10 34 ergs -1 for the 
Crab Nebula (assuming a distance of 2 kpc). Thus, the effi- 
ciency of VHE 7-ray production relative to the spin-down 
rate (L/E) of PSRB1509-58 (~0.6%) is much larger than 
that of the Crab pulsar (~0.016%). Such high efficiency 
has been proposed previously IjAharonian et al. 1997|) . 

A simple one- zone IC model (Khclifi 2002) can be 
used to reproduce the VHE spectrum of the whole neb- 
ula. Therein, a population of accelerated electrons with 
a power-law energy spectrum is assumed. These particles 
lose their energy by synchroton emission in a magnetic 
field, assumed to be uniform within the PWN, and by IC 
scattering on seed photons. The electron spectrum and 
the magnetic field strength are adjusted such that the 
computed synchrotron and IC spectra match the X-ray 
and VHE 7-ray observations. The energy spectrum of the 
whole nebula measured by BeppoSAX (|Mineo et al. 200l)l 
is used for this fit. By using only the cosmic microwave 
background radiation (CMBR) as seed photons, it is not 
possible to reproduce the VHE spectrum. Thus, infrared 
(IR) photons from dust and starlight are added, using a 
recent parametrisation of the interstellar radiation field 
(ISRF) used for the GALPROP code UStrong et al.~2000| ). 
The energy density of the dust and starlight components 



of the ISRF are however kept free to account for possible 
local variations. 

The best fit model is represented in Figure in the 
VHE band with the contribution of the different seed pho- 
tons to the IC spectrum. The fitted spectral index of ra- 
diating electrons is 2.9. The fitted mean magnetic field 
strength is ~17 /iG, which is about a factor of two higher 
than the lower limit estimated previously (Gaensler et al. 
2002 ) . The best fit energy density of the dust component is 
about 2.3eVcm -3 which is more than twice the nominal 
value of the GALPROP ISRF. This apparent enhance- 
ment of dust density may be due to the material swept up 
by the SNR G 320.4-1.2 (du Plessis et nl. HMfy . As the en- 
ergy density of the starlight component is unconstrained, 
it is fixed to the GALPROP value (~1.4eVcm -3 ). 

Although this spectral model is compatible with all 
available measurements, the derived value of B, when 
combined with the pulsar spin-down age of 1700 years, 
suggests the existence of a cooling break in the electron 
spectrum of ~24TeV. This corresponds to a photon en- 
ergy above ~3 TeV in the IC spectrum, depending on the 
energy of the seed photons, which should be detectable. 
Nonetheless, if the age of the system is closer to 20 000 
years as suggested by Gvaramadze l)2001[) . the above spec- 
tral model would be self-consistent. It should also be noted 
that there is noticeable correlation between the fitted B 
value and the photon energy density from dust. Thus, if 
this local density was higher, the corresponding B value 
would increase and thus the cooling break energy in the 
IC spectrum would decrease. 

In summary, this detection represents the first morpho- 
logical identification of an extended pulsar wind nebula in 
the VHE 7-ray band. The H.E.S.S. observations provide 
more direct information on accelerated particles than soft 
X-ray data, due to the uncertainties in the estimation of 
the magnetic field and its spatial variations. 
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